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The Casimir effect in a dispersive and absorbing multi- 
layered system is considered adopting the (net) vacuum-field 
pressure point of view to the Casimir force. Using the prop- 
erties of the macroscopic field operators appropriate for ab- 
sorbing systems and a convenient compact form of the Green 
function for a multilayer, a straightforward and transparent 
derivation of the Casimir force in a lossless layer of an other- 
wise absorbing multilayer is presented. The resulting expres- 
sion in terms of the reflection coefficients of the surround- 
ing stacks of layers is of the same form as that obtained by 
Zhou and Spruch for a purely dispersive multilayer using the 
(surface) mode summation method [Phys. Rev. A 52, 297 
(1995)]. Owing to the recursion relations which the general- 
ized Fresnel coefficients satisfy, this result can be applied to 
more complex systems with planar symmetry. This is illus- 
trated by calculating the Casimir force on a dielectric (metal- 
lic) slab in a planar cavity with realistic mirrors. Also, a 
relationship between the Casimir force and energy in two dif- 
ferent layers is established. 

PACS numbers: 12.20.Ds, 42.60. Da 



I. INTRODUCTION 

Originally, the Casimir effect was predicted as a feature 
of the electromagnetic field between two neutral ideally 
conducting plates and consisted in the appearance of an 
attractive force between the plates. The force is due to 
the change of the zero-point energy of the field in the con- 
fined space In this special case, however, the Casimir 
force can also be viewed as the long-range van der Waals 
force It becomes appreciable in the submicron 

and rapidly increases in the nanometer range. As such, it 
may strongly affect processing in nanotechnology as well 
as functioning of micro- and nanomachines and devices 
1^. Clearly, these new developments pose the problem 
of realistic calculations of the Casimir force on objects in 
complex enviroments. 

In contrast to the highly idealized system considered 
by Casimir Lifshitz Q calculated the force between 
two thick (semi-infinite) dielectric slabs by taking into ac- 
count the dispersion and absorption in the dielectrics as 
well as the temperature effects. In this respect, his the- 
ory is far more realistic and, as the effects of finite con- 
ductivity and dissipation in the metal can be observed 
in the recent high-precision experiments his re- 

sult for the force at zero temperature is standardly used 
when analyzing the Casimir force in the planar geometry 
|^,|o|. The Lifshitz approach is based on the calculation 



of the electromagnetic field due to the randomly fluctuat- 
ing currents in the dielectric slabs and on the subsequent 
calculation of the Maxwell stress tensor in the region in- 
between. Owing to its complexity, however, it has never 
been extended to the calculation of the force between 
multilayered stacks although the generalization of the fi- 
nal result to this configuration is fairly obvious. 

The Casimir effect in multilayered systems is usu- 
ally considered using either the surface mode summation 
method 0Jl|,|l|,|l|l (see also Refs. |,|) to calculate the 
change in the electromagnetic field zero-point energy due 
to the presence of the dielectric stacks, or the stress ten- 
sor method to calculate directly the vacuum field 
pressure on the stacks. Strictly speaking, the mode sum- 
mation method applies only to purely dispersive (lossless) 
systems as only in this case the mode frequencies are real. 
However, when expressed as an integral over the imag- 
inary frequency, the final result for the Casimir energy 
(force) turns out to be applicable to absorbing systems 
as well. An indication that this must be so is the fact 
that the dielectric function is always real on the imagi- 
nary axis irrespective of whether the system is absorbing 
or not Thus, while in their calculation of the force 
in a multilayer Zhou and Spruch [p^ assumed a purely 
dispersive system, Klimchetskaya et al. |Q recently con- 
sidered a similiar but absorbing system. 

On the other hand, being a local approach, the stress 
tensor method does not necessarily imply a lossless sys- 
tem. Since the stress tensor cannot be defined macro- 
scopically for an absorbing medium the only nec- 
essary assumption is actually that the region where 
the vacuum-field pressure is calculated is nonabsorbing, 
whereas the other parts of the system may generally 
be dissipative. Despite this fact, numerous papers in 
the past used the stress tensor method to calculate the 
Casimir force assuming, at most, a dispersive but nonab- 
sorbing system. One of the reasons for that is certainly 
lack of knowledge on the proper form and properties of 
macroscopic field operators aproppriate for an absorbing 
system at that time. 

The first calculation of the Casimir force between two 
absorbing slabs is due to Kupiszewska who mod- 
eled dielectric atoms as a collection of harmonic os- 
cillators coupled to a heat bath that absorbs energy. 
Only the modes propagating normally to the slabs were 
considered, so that this approach was effectively one- 
dimensional (ID). Describing the reservoir through a 
damping constant and the Langevin force and solving for 
the field operators, Kupiszewska obtained for the force 
between the slabs the same expression in terms of their 



1 



reflection coefficients as that obtained previously for an 
inert jl^] or a lossless Q ID system, except that this 
time the dielectric function of the slabs was complex. Re- 
cently, this result was rederived using a Green function 
method for quantizing the macroscopic field in (ID) ab- 
sorbing systems in conjuction with a scattering matrix 
approach |2|] and was also extended to two identical ab- 
sorbing superlattices [ p2| . Very recently, Esquivel-Sirvent 
et al. demonstrated an alternative Green function ap- 
proach that makes the quantization of the field within 
the slabs unnecessary and calculated the Casimir force in 
an asymmetric configuration |23|| which was earlier con- 



sidered only in the lossless case |2C 

Owing to their complex structure, an explicit calcu- 
lation of the field operators as in Refs. | l^ , |l9| , po| , ^ 
is highly impractical in the general case of a three- 
dimensional (3D) dissipative inhomogeneous system. 
However, as pointed recently by Matloob (see also 
Rcf. [pT| ), using the fluctuation-dissipation theorem and 
the linear response theory, the field correlation functions 
needed to calculate the stress tensor can be expressed in 
terms of the (classical) Green function for the system. In 
this way, only the knowledge of the Green function is, 
therefore, actually needed to calculate the Casimir force. 
Using this method, Matloob and Falinejad recently con- 
sidered the Casimir force between two identical absorbing 
dielectric slabs |2^] . Very recently Mochan et al. \ 26 1 gen- 
eralized their Green function method to three dimen- 
sions and calculated the Casimir force between two arbi- 
trary slabs. Expressing the refiection coefficients of the 
slabs through the generalized surface impedances, these 
authors argued that their formal result could be applied 
to rather general but not chiral media, also including non- 
local inhomogeneous dissipative slabs. In Refs. [ p5|j26| 
the space between the slabs was assumed empty. 

In this work we calculate the Casimir force in a lossless 
dispersive layer of an otherwise absorbing multilayer by 
employing the macroscopic field operators as emerge from 
a recently developed scheme for quantizing the electro- 
magnetic field in inhomogeneous dissipative 3D-systems 
|p7| , p8| and using a convenient Green function for a mul- 
tilayer 1^ . In this way, we obtain a general result for the 
Casimir force in stratified local media. In addition, us- 
ing the properties of the generalized Fresnel coefficients, 
we derive a relationship between the Casimir force and 
energy in two different layers and demonstrate the ap- 
plicability of the theory to more complex planar systems 
by calculating the Casimir force on a diectric slab in a 
realistic planar cavity. 



II. THEORY 

Consider a multilayered system described by the di- 
electric function e(r,w) = e'(r,a') -I- ie"(r,w) defined in 
a stepwise fashion, as depicted in Fig. 1. The Casimir 
force in a layer corresponds to the net vacuum-field pres- 



sure in the multilayer with respect to the pressure in the 
infinite layer (medium). Accordingly, the force 'Pj/i on 
a stack o f lay ers that separates a jth and an /th layer is 
given by [|o| 



(2.1) 



where Tj^^z is the zz component of the regularized stress 
tensor in the jth layer 



^j,ab — -'-j.ab -'^j,abi 



(2.2) 



with Tj^ab and ^f, being the corresponding Maxwell 
stress tensors in the multi laye r and in the infinite medium 
(j), respectively. In Eq. ( ^.l[ ), A is the area of the stack 
and the -I- (— ) sign apphes if Z > j (^ < j). Since the 
regularized stress tensor vanishes in the outmost layers, 
we have for = fj/o and fj+ = /j/„: 



(2.3) 



so that Tj^zz coincides with the force per unit area acting 
on the left (right) stack of layers bounding the layer (j). 
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FIG. 1. System considered schematically. The dashed line 
represents the plane where the stress tensor is calculated. 

Replacing fleld variables in the classical Maxwell stress 
tensor |l^] by the corresponding Heisenberg operators 
and taking its average |l^], Tj ^z in a lossless layer (j) 
[e"(tj) = 0] is given by 



{EzDz 



E||.D|| 



BzHz — B|| • H||) 



re(i) ' 
(2.4) 



where we have suppressed the argument (r, t) of the field 
operators and the brackets denote the expectation value 
in the vacuum state of the field. In order to calculate 
the correlation functions that appear in Eq. ( 2.4) 



we 



use the properties of the macroscopic field operators ap- 
propriate for absorbing systems [ p7| . These operators 
are decomposed into their "annihilation" and "creation" 
components according to 
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E(r,t)=/ duj'E{r,uj)e-"^* + H.c. (2.5) 

and, with the constitutive relations 

D(r, oj) = £(r, tj)E(r, uj) + 47rPjv(r, uj), 

B(r,c^) =H(r,c^), (2.6) 

obey the standard macroscopic Maxwell equations. Here 
Pjv(r,Ci;) and 'P\^{r,uj) are the noise polarization opera- 
tors related to the dissipation in the system and obeying 
the commutation rules (in the dyadic form): 

[P^(r, u;), Pjv(r', c')] - ^'"}Y^ IS{r - v')6{cu - u;'), 



47r2 



(2.7) 



where I is the unit dyadic. Therefore, any (annihila- 
tion) field operator is related to PAr(r, uj) via the classical 

Green function G(r, r';a;) satisfying 



V X V X -e{r.uj)—I- 



G(r,r';w) = 47rI5(r- r') 



(2.8) 



according to 
E(r,w) = 



d^r'G(r,r';w) •PAr(r',w). (2.9) 



As a consequence, all field correlation functions can be 
expressed through the Green function in accordance with 
the fluctuation-dissipation theorem |31|. In particular, 
for the electric-field correlation function we have [ETI 



(E(r,L^)Et(r',L^')) = ^ImG(r, r'; w)J(c^ - c^'), 



and the magnetic-field correlation function is easily ob- 
tained from this expression using B(r, w) = (— ic/a;)V x 
E{r,uj). ^ 

Applying the above results to the jth layer and taking 
into account that ej{uj) is real and that PAr(r,w) =0 in 
this region, we find for the relevant correlation functions 
in Eq. (0|): 



(E(r,t)D(r,i)) 



re(j) 



dijjfc|(tj)ImGj(r, r; lo), 



(2.11a) 



(B(r,t)H(r,t)>re(^.) --y^ dculmGf {r,r;uj), (2.1 



lb) 



where kj (uj) = yj£j{uj)ijj /c is the wave vector in the layer, 

Gj(r,r';u;) is the Green function element for r and r' 
both in the layer (j), and 



Gf(r,r';tj) = V X Gj(r,r';tj)x V' 



(2.12) 



is the corresponding Green function element for the mag- 
netic field . With the above equations inserted in Eq. 
(p. 4), the stress tensor Tj^^z is expressed entirely in terms 
of the Green function and its derivatives analo gou sly to 
Eq. ( P^sl) below. Similarly, through Eqs. ( ^.4| ) and 
( ^.ll| ) applied to the infinite medium (j), the stress ten- 
sor T?^^ is given by the same expression with the infinite- 
medium Green function G"(r, r';a;). Therefore, the reg- 
ularized stress tensor Tj^^z is expressed as 

fi 



T- — 

J^3,zz — 



47r 
G 



-Im 



duo 
2 



(r,r;w) - qyr,r;w) 



^(r,r;^)-Gff (r,r;a;)}, (2.13) 



where 



Gf (r, r'; oj) = G, (r, r'; cj) - G°(r, r'; lj) (2.14) 

is the Green function for the scattered field in the jth 
layer and G^^^ (r, r'; lu) = Gf^^{r, r';LLi) + Gfyy{r, r'; lu) is 
its parallel trace. 

A convenient form of G^'^(r,r';w) for a general multi- 
layer is derived in Ref. |29]. In the Appendix we quote 
this Green function and calculate the expression in the 
curly brackets of Eq. ( 2.13| ). Inserting Eq. (|A8|), wc find 



T- — — - 



■Re 



duj 



(27r)2 
dkkn 



q=p,s 

'^~Dqj{i^,k) 

q=p,s 



Dqj{uj,k) 

"^■^"^'■"^ (2.15) 



Dqj{ii,k) 



(2-10) where /3j(w, k) = ^~^{uj)-k^, 



Dqj (tj, fc) = 1 - (w, fc)r^ , (u), k)e 



2iPjdj 



(2.16) 



and rj^(u;,fc) are the reflection coefficients of the right 
and left stack of layer s bounding the jth layer. The sec- 
ond line in Eq. ( 2.15| ) has been obtained by converting 
the integral over the real w-axis to one along the imagi- 
nary a;-axis in the usual way, letting uj ~ i^. 



k) = k) = i^e,izOe/c^ + k\ (2.17) 

and noting that the expression in the brackets is real on 
the imaginary axis. We see that the regularized stress 
tensor is uniform across the layer. Although expected 
on invariance grounds, this is not a trivial result and, as 
is clear from the derivation in the Appendix, it is due 
to cancellation of the z-dependent terms in the electric 
and magnetic contributions to Tj,zz irrespective of the 
dielectric properties of the surrounding stacks. 
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Knowing the force, the Casimir energy £j in the layer 
can be calculated using 



with the condition that £^ 



for di 



( p. 151) and (gj), we find 



(2.18) 
oo. From Eqs. 



£j — him 



i2n) 



OO 



(27r)2 

OO 

dkk 



E 

q=p,s 

E 

q=p,s 



In Dgj{uj,k) 



\nD,j{t^,k). (2.19) 



This e qua tion, a s wel l as that for the force [combined 
Eqs. (2.3) and ( 2.15| )], agrees in form with the corre- 
sponding result of Zhou and Spruch fl^ derived using 
the (surface) mode summation method and starting from 
a simple model of a purely dispersive multilayer. How- 
ever, in this work the contributions of all (propagating 
and evanescent) modes are naturally taken into account 
on an equal footing through the Green function. Fur- 
thermore, since the Green function employed refers to a 
general absorbing multilayer, so do the obtained results 
except, of course, for the region where the Casimir force 
is calculated. 

The Casimir energy and force vary in a stepwise fashion 
across the multilayer and we end this section by pointing 
out a relationship that exists between their values in two 
different layers, say, layers (j) and (l). Indeed, assuming 
that I > j, for example, and using recursi on re lations for 
the reflection coefficients given by Eq. ( A3a ), one may 
prove that the followin g r elation exists between the D 
functions for the layers |29]: 



Dqi{l - r:„r:^e 



1 ^9 p2i/3jdj 



DqAl 



^1 „<1 „2ifi,di 



(2.20) 



Combining this with Eq. ( 2.19 ), we find that the respec- 
tive Casimir energies are related according to 

fi 



£i = £, 



q=p,i 



(2^)2 
1 

In ■ 



/ dkk X 
Jo 



1 



rj^j(a;, fc)rj_(aj, k)e'^'^^i'^i 



(2.21) 



A similar relation is obtained for the forces in two layers, 
but the resulting expression is not particularly illuminat- 
ing unless Ej = Such a situation arises, for example, 
when a planar object is embedded in a planar cavity. In 
this case, we find 



n 



d^ 



dkkn 



1 



Dqjii^,k) 



^-rj/iii^,k)r]_{uj,k)e 



2if3dj 



1 - 



r«/^.(w,fc)rf^(tj,fc)e2*/3di 



ijj—i^ 



where P (k) is the perpendicular wave vector in both 
layers. 



III. DISCUSSION 

Most of the previously obtained results for the Casimir 
force and energy in a specific planar configuration are re- 
covered from the results derived in the preceding sec- 
tion simply by specifying the corresponding reflection 
coefficients and material parameters. Thus, for exam- 
ple, the results for the three-layer (£1,63, £2) configu- 
ration considered by Lifshitz |^1[ are obtained letting 
£j = £3, r^i and rj_(^ — > r^j, where r'^- are single- 

interface refiection coefficients given by Eq. (A4), and 
the results for the five-layer (£4, £1, £3, £2, £5) configura- 
tion considered by Zhou and Spruch ||l3|] are obtained 
letting £j = £3, — > r|j^4 and rj_|_ ''325; where the 
three-layer re flect ion coefficients are obtained from recur- 
rencies Eq. ( |A3a ). Similarly, the results for the system 
consisting of two identical slabs, recently considered by 
Matloob and Falinejad (251, are obtained letting £, = 1 



and r' 



3± 



f', where r'^ are the reflection coefficients of 



a symmetrically bounded slab [see Eq. (3.5) below], etc. 
Specially, the Casimir force and energy in a (dispersive) 
planar cavity formed by two ideally reflecting (conduct- 
ing) slabs are obtained with i"j_r'j^ = 1. 

We also note that these equations correctly reproduce 
the corresponding results which emerge from the ID con- 
siderations. Indeed, taking only the fc = contribution 



in Eq. (2.15), we flnd from the flrst line in that equation, 
for example. 



'd.k.ic.) ^ (3.1) 



where Dj{uj) = Dqj{uj,0) [Eq. ( |2l^ )] is the same for 
both polarizations. With a simple algebra, this equation 
can be rewritten as 



rplD _ ^ 
j.,zz 



dujkj{ijj) X 



|1 - rj_(w)rj+(tj)e2«'='3('^)'^j|2 



(3.2) 



which is in accordance with the Casimir force obtained by 
several authors for the respective systems they considered 
]ll,||j20|||,||,||. 

Owing to the recursion relations which the generalized 
Fresnel coefficients satisfy, the obtained results can be 
applied to more complex systems with planar symmetry. 
As an application of the theory, we illustrate this by de- 
riving the Casimir force on a dielectric, or a metallic, slab 
[dielectric function Es, thickness I] in a cavity [dielectric 
function e] with realistic mirrors [reflection coefficients 
r\ and rf], as depicted in Fig. 2. The force on the slab 
/ = /2_— /i_ in this configuration can be calculated from 



(2.22) Eq. ( |2.22| ). The function Dqi [Eq. ( |2.16| )] is straightfor 



wardly obtained letting r\_ 



- ^1 



and 



"2+ 



and using 



Eq. (A3a) to determine the reflection coefficients r\_^. 
We find (the polarization index q is omitted) 
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D^-^-^^ir+J^^^^)e^'''^. (3.3) 

where r = ri/2 = ''2/1 Etnd i = ti/2 = t2/i are Fresnel 
coefficients for the slab. 




L/2 



da 



Mo 



FIG. 2. A dielectric slab in a planar cavity shown schemat- 
ically. The arrow indicates the direction of the force on the 
slab. 

This gives 



f = — - / df dkkn x 



E 

q=p,s 



N 



uj—i^ 



(3.4) 



where the expression in the brackets is to be calculated 
for g-polarization. Using Eqs. (A3) and (|A4|), r and t 
can be further expressed entirely in terms of the reflection 
coefficient for the cavity-slab interface p = {1 — rj) / {1 + rj) 
[where rjP = e(3s/£sf3 and 77* — f3s/P] as 



1 



t = 



(1-pV 

I _ „2^2if3J 



Note that for a perfectly conducting (sg 
have rP = -r" ^ 1 and i« = 0. 



(3.5) 

00) slab, we 



The force /, as given by Eq. (3.4), may be positive 
or negative, depending on the dielectric properties of the 
slab and cavity mirrors as well as on the position of the 
slab. One may easily verify that this equation gives the 
correct result for the force on a perfectly conducting plate 
in an empty cavity with ideally reflecting walls. Indeed, 
since in this case P — Q,N factorizes and Eq. (3.4) splits 
to (r^ = rri = rr2 = rir2 = 1) 



f = f-id2) - f-{di), 



(3.6) 



where 



37r^c^ 



dec 



1 



(3.7) 



is the force on the left mirror of a dispersive ideal cavity 
[cf. Eqs. (2.3) and (2.15), with the index j dropped]. 
The second line here is obtained upon a partial inte- 
gration over C and upon calculating the ^-derivative of 
the integral over k (see Rcf. [ p2[). F or the empty cavity 
[s{i^) = 1], the integrals in Eq. (3/7) become elementary 
giving the well-known result 



/ 



240 



1 



(3.8) 



according to which the plate is attracted to the closer 
cavity mirror. For a partially transmitting plate, the 
vacuum-field fluctuations in regions (1) and (2) of the 
cavity are no longer independent of each other and con- 
siderable deviations from the above simple result may 
occur especially for realistic cavity mirrors. Clearly, in 
this case, in order to explore the combined effect of the 
nearby walls on a planar (nano)object, one must analyze 
Eq. (3.4) numerically. 



IV. SUMMARY 

Using the properties of the macroscopic field opera- 
tors appropriate for dissipative systems and a convenient 
Green function for a multilayer, in this work we have ob- 
tained general results for the Casimir force and energy 
applicable to local layered absorbing systems. We have 
also established a relationship between the Casimir force 
(and energy) in two different layers and, as an application 
of the theory, calculated the Casimir force on a dielectric 
slab in a realistic planar cavity. 
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APPENDIX: GREEN FUNCTION 

Denoting the (conserved) wave vector parallel to the 
system surfaces by k = (fc^;, fcj,), we write the wave vector 
of an rightward (leftward) propagating wave in an /th 
layer as Kt*^ = k ± /3;z, where 



A = V^i'-fc'^A' + ^A", A'>o, A">o. (ai) 
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With this notation, the Green function dyadic for the 
scattered field in the jth layer reads |29 



+,J^e^ft^+e^-^.(k) [e-^.(_k)e-''^-' +rJ_e+.(-k)e^'^^-' ] } 

z_ = z, z+ = — z, < 2 < (ij, 

ej.(k) = ^(±/?,k + fcz), e5(k)=kx£ = n, (A2) 

where r|_|_ = ''j/„(o) ^^re, respectively, the transmission 
and reflection coefficient of the upper (lower) stack of 
layers bounding the layer (j). Clearly, for the outmost 
layers, ^ = n (0), we have r^^^ = and = 0. Also, 
one must let dn {do) = since these quantities appear 
only formally. The remaining Fresnel coefficients satisfy 



q ^ + VJ j/' j/fc 

'i/j/k i/i 1 _ ^9 g2i0jdj 



(A3a) 



and G^- '^''^(r, r; w) is given by this equation multiplied 
by and with p ^ s. The traces G^'^n (r, r; w) and 
G^l'l'*'^(r,r;a;) can be easily recognized from these equa- 
tions and one has, for example. 



G-||(r,r;c.)-G-,(r,r;a; 

0' 



27r J fij 



D 



-,2if3j zj. 



while G^|p(r, r; w) — G^'^^(r, r; w) is given by this equa- 
tion with p s. Adding these two quantities, one finds 
that the curly bracket in Eq. ( 2.13 ) is equal to 



{•■•} = 



d^k 



(3, X 



(2^)2^^ 

„P JP „2i0jdj s s „2if3jdj 

1- I ' j-' 3 + ^ 



(A8) 



'i/s/k i_^q^^^q^^^2ii3,d, fc/jA' 
and, for a single i — j interface, reduce to 



_ |3^ - lljP, , 



ij r^3 ,q I q / -i i q \ 



(A3b) 



(A4) 



where 7? = £i/£j and 7^"^- — 1, respectively. 

Performing the derivations indicated in Eq. ( 2.12 ) and 
using 

K±(k)xe±(k)=fc,C,e±^.(k), p' = s, s' = p, (A5) 

we find that Gf'^'^{r,r';u!) is given by Eq. (A2) mul- 
tiplied by — fc| and with e^- ^f'j- Noting that the 
equal-point Green function dyadics consist only of diag- 
onal elements, we easily find 



Gr(r,r;^) 



d^k 



(A6) 



^P3 } 
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